Effect of tensile stress on the in-plane resistivity anisotropy in BaFe 2 As 2 
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The effect of uniaxial tensile stress and the resultant strain on the structural/magnetic transi- 
tion in the parent compound of the iron arsenide superconductor, BaFe2As2, is characterized by 
temperature-dependent electrical resistivity, x-ray diffraction and quantitative polarized light imag- 
ing. We show that strain induces a measurable uniaxial structural distortion above the first-order 
magnetic transition and significantly smears the structural transition. This response is different 
from that found in another parent compound, SrFe2As2, where the coupled structural and mag- 
netic transitions are strongly first order. This difference in the structural responses explains the 
in-plain resistivity anisotropy above the transition in BaFe2As2. This conclusion is supported by 
the Ginzburg-Landau - type phenomenological model for the effect of the uniaxial strain on the 
resistivity anisotropy. 



INTRODUCTION 

At ambient conditions, the parent compounds of iron- 
arsenide superconductors, AFe2As2 (A = Ba, Ca or Sr), 
crystallize in the tetragonal ThCr2Si2 structure [TJ [2]. 
On cooling, they undergo a structural phase transition 
with the lattice symmetry lowered from tetragonal to or- 
thorhombic at a characteristic temperature Tto (170 K 
for A=Ca [J, 0, 205 K for A=Sr [3j and 135 K for 
A=Ba PQ. We denote the compounds as A122 in the fol- 
lowing). This transition is accompanied or followed by 
long-range magnetic ordering into an antiferromagnetic 
(AFM) stripe phase at the Neel temperature, T/v 0]. In- 
deed Tat = Tto in compounds with A=Ca [5 and A=Sr 
[6], where the transition is sharp and strongly first order. 
In BaFe2As2 Tjy < Tto-, and the structural transition is 
second order, whereas the AFM transition is first order 

Since the doping (or pressure) - dependent supercon- 
ductivity in A122 iron arsenides exhibits the highest T c 
close to the point of complete suppression of the struc- 
tural/magnetic order, understanding the mechanism of 
these transitions is very important for understanding the 
nature of superconductivity. The parent compounds of 
iron arsenides are metals, so it is suggested that their 
magnetism is of itinerant character due to a spin density 
wave (SDW) instability of the multi-band Fermi surface 
[H[T0] • On the other hand, it has also been suggested that 
the magnetism can arise in a local moment picture due to 
magnetic frustration [11] and/or orbital ordering [T2HT5] . 
Therefore, it is important to conduct measurements that 
characterize the normal state anisotropy of the electronic 
structure in the vicinity of T/v- From first principles cal- 
culations, the electronic anisotropy of iron pnictides was 
predicted to be fairly high in the orthorhombic ab— plane 
below T/v [T6H20] . On the other hand, ARPES measure- 
ments suggest that a notable energy splitting between d xz 
and d yz orbitals appears below the transition [2TJ [22] . 



An insight into intrinsic anisotropy became possible 
after the development of detwinning techniques, using 
uniaxial tensile p~6j [23] or compressive stress [24] [25] 
(see Ref. 26 for a review). Electrical resistivity measure- 
ments in the detwinned state found the in-plane resis- 
tivity anisotropy to have an unusual temperature depen- 
dence of the ratio Pb/p a , peaking just below Tto with 
maximum pb/p a = 1.2,1.4,1.5 for A=Ca, Sr, Ba, re- 
spectively [16] [23]. Surprisingly, in BaFe2As2 the high 
temperature "tail" of the anisotropy is found even in the 
nominally tetragonal phase above T T o- The anisotropy 
above Tto becomes most pronounced in slightly doped 
Ba(Fei_ a ,Co :c )2As2 (BaCol22), x < 0.06[24], but is not 
observed in either Cal22 or Srl22. This difference be- 
tween Cal22 and Srl22 compounds on one hand, and 
Bal22 and BaCol22 compounds on the other, may be 
related to the sharpness of the structural phase transi- 
tion, i.e., a strongly coupled first order [5], [6] vs. split 
second order transitions [7] [8] . 



In this paper, we study the evolution of the electronic 
and structural anisotropy of detwinned Bal22 with spe- 
cial attention to the effects of the applied strain required 
to detwin the samples. With a much better strain con- 
trol, we find that the effect above the transition arises 
from an anomalously large structural susceptibility of the 
crystals to the applied strain. This strain further sepa- 
rates the already split structural and magnetic transi- 
tions in BaFe2As2, as found in the detailed analysis of 
the temperature-dependent x-ray and polarized optical 
imaging. We model the effect of the applied strain by us- 
ing a phenomenological Ginzburg-Landau - type model 
and show that the difference in the response is directly 
linked to the order of the magnetic/structural transitions. 




FIG. 1. (Color Online) The top right panel shows a schematic 
of the horseshoe with the potential leads used to apply the 
strain by adjusting the push-screw. The current leads are 
mounted strain-free. The top left panel shows the photograph 
of an actual sample with soldered contacts. The area "A" on 
the left side of the sample represents the unstrained region. Its 
polarized light image at 5 K (bottom left panel) reveals clear 
domain pattern by the alternating blue and orange stripes. 
The region "B" is located between the potential contacts in 
the strained part of the sample. It is shown in its detwinned 
state in the bottom right panel. 

EXPERIMENTAL METHODS 

Single crystals of Bal22 were grown from an FeAs flux 
as described in Ref. [27| Crystals were cut into strips with 
typical dimensions of 0.6 mm wide, 3 to 5 mm long and 
had a thickness of approximately 0.1 mm. The lengthwise 
cuts were made parallel to the tetragonal [110] direction, 
which will become the orthorhombic a Q or b Q axes be- 
low the transition temperature. Cutting directions were 
estimated by eye using polarized optical images of the 
domain structure and natural facets on the sample. 

Polarized light images were taken at temperatures 
down to 5 K using a Leica DMLM polarization micro- 
scope equipped with a flow- type 4 He cryostat, as de- 
scribed in detail in Ref. [28j High-resolution images were 
recorded with a spatial resolution of about 1 jam. Mea- 
surements were done with the polarizer and analyzer 
nearly crossed. 

The leads for electrical resistivity measurements were 
formed by Ag wires, 50 fim in diameter, soldered to 
the samples with tin [29 . A photograph of the sam- 
ple with wires is shown in the top left panel in Fig. [I] 
Four - probe measurements were conducted in a Quan- 
tum Design PPMS from 5 K to 300 K. Measurements 
were first made on a free standing sample, and then the 
voltage leads were attached to a horseshoe, as schemat- 
ically shown in the top-right panel of Fig. [T] while the 
current leads were mounted so as to produce no strain 
on the sample. The strain was applied by means of a 




FIG. 2. (Color Online) Left panels show polarized light im- 
ages of the strained portion "B" of the sample below (130 K, 
top) and above (160 K, bottom) the structural transition tem- 
perature. The intensity in red, green and blue (RGB) chan- 
nels for each pixel was digitized using 256 intensity bins. The 
total intensity was found by summing the intensities of all 
pixels in the selected region. The RGB color was character- 
ized by the percent contribution of each channel to the total 
intensity, plotted against the value of each bin to produce the 
histograms in Fig. [2] Right panels show the RGB histograms 
of a small area of the strained portion "B" of the sample, in- 
dicated by the red square in the top left panel. Whereas the 
blue channel remains almost unchanged, the intensities of the 
green and red channels shift dramatically indicating overall 
spectral change. The temperature dependence of this change 
is quantitatively analyzed in Fig. [3] 

stainless push screw expanding the legs of the horseshoe. 
The temperature-dependent resistivity was measured af- 
ter every strain increment. For the evaluation of the ten- 
sile stress value we compared our data with the data of 
T. Liang et al [30], who found a roughly +5 K shift of 
the tetragonal-to-orthorhombic transition feature upon 
a stress change from 15 to 50 MPa. In our case a total 
shift of approximately 3 K was achieved in four equivalent 
stress increments, which suggests that each strain incre- 
ment is in the 4-5 MPa range for our horseshoe straining 
device. Therefore the strain at the highest level is esti- 
mated to be in the range of 16-20 MPa. 

The sample was periodically imaged via polarized mi- 
croscopy. The bottom panels in Fig. [T] show two areas 
of the sample: area A (left panel) is located between 
current and potential leads and remains twinned during 
measurements; area B (right panel) is located between 
the potential contacts in the strained part of the sample 
and becomes nearly twin-free under strain. The appli- 
cation of uniaxial stress makes it energetically favorable 
to align domains with their long a Q axis along the strain, 
giving rise to an increasing volume fraction of one domain 
orientation above the rest. 
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FIG. 3. (Color Online) Normalized temperature variation 
of the green color channel's intensity through the magnetic 
transition in detwinned crystals of BaFe2As2 and CaFe2As2. 
Arrows indicate the temperatures of the respective magnetic 
transitions. While there is a clear signature of the transition 
at about 135 K in BaFe2As2, the curve changes smoothly 
through the transition with a second order character. In 
CaFe2As2, the transition is quite sharp around 160 K, and 
is strongly first order. 

The detwinned crystals were studied by high-energy 
x-ray diffraction in the MU-CAT sector (beamline 6ID- 
6) at the Advanced Photon Source at Argonne National 
Laboratory. Measurements were made with X-rays at 
99.3 keV, giving an absorption length of approximately 
1.5 mm and thus allowing full penetration of the typi- 
cally 0.1 mm thick samples. Entire reciprocal planes were 
recorded with the c-axis parallel to the incident beam. 
During measurements, the samples were rocked through 
two independent angles perpendicular to the beam (see 
Ref. 31). The direct beam was blocked by a beam stop 
and diffraction images were recorded by a MAR345 im- 
age plate positioned 1680 mm behind the sample. The 
beam size was reduced to 0.2x0. 2mm 2 by a slit system. 



RESULTS 
Polarized Microscopy 

The unit cell in BaFe2As2 doubles in size and rotates 
by 45 degrees upon cooling through the tetragonal to 
orthorhombic transition. This leads to the formation of 
the domain walls at 45 degrees with respect to the sample 
edges (see bottom left panel of Fig.[T]). The orthorhombic 
a Q and b Q axes inside the domains are at 45 degrees to 
the twin boundaries (see Ref. [28]). Therefore the highest 
contrast of domain imaging is achieved when the sample 
is aligned with a long [110] tetragonal direction at 45° to 
the polarization direction of the linearly polarized light 
(parallel and perpendicular to the orthorhombic a Q in 
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FIG. 4. (Color Online) Temperature evolution of the two- 
dimensional x-ray diffraction pattern near the tetragonal 
(220) Bragg diffraction peak. Left and right columns of im- 
ages show diffraction patterns in the unstrained and strained 
parts of the crystal, respectively. Four spots in the unstrained 
part at 6 K (top left) represent four domains in the sample 
with domain populations (proportional to integrated inten- 
sity) ranging between 19 and 31%, close to random (see also 
Ref. 1 28 p . In the strained portion of the sample (6 K, top right 
panel), the dominant domain occupies nearly 90 percent of 
the volume of the sample area probed by the x-ray beam. Be- 
tween 132 and 134 K, the sample undergoes an orthorhombic 
to tetragonal structural transition. The second-order nature 
of the transition is evidenced by the lack of coexistence of 
orthorhombic and tetragonal domains. This coexistence is 
clearly observed in Srl22 (see Ref. [23]), in which the tran- 
sition is strongly first-order, see Fig. [5] below for schematic 
elaboration. 

different domains.) 

The optical contrast of the domains is determined by 
the anisotropy of the bi-refiectance and is proportional to 
the orthorhombic distortion. It increases with decreasing 
temperature, and is weaker in Bal22 than in Cal22 [28] . 
Simultaneously, due to dispersion of the bi-reflectance, 
initially white incident light on reflection acquires color 
depending on the orientation of the orthorhombic axes 



4 




(040) ( 



FIG. 5. (Color Online) Schematic diagrams of the displace- 
ments of atoms in the twinned orthorhombic phase and the 
resulting Bragg reflections. As demonstrated in the top panel, 
a perfect crystal with equal populations of each domain ori- 
entation results in a square pattern between the (400) and 
(040) orthorhombic reflections. Conversely, the bottom panel 
illustrates the result of an unequal distribution of domain ori- 
entations. Here the angle between the 04 and the 02 domain 
orientations is significantly smaller and consequently moves 
the reflections closer together. Further, the population of each 
domain is proportional to the intensity of its Bragg reflection. 
These effects can be seen in the X-ray data of Fig. [4] especially 
the T=132 K panels. 



with respect to the polarization direction of the incident 
light. This results in different colors of structural do- 
mains as can be seen in the bottom left panel of Fig. [T] 
Therefore, the color of the image contains information 
about the orthorhombic distortion, albeit in arbitrary 
units, and can be used for the analysis of its tempera- 
ture dependence even in the detwinned state below the 
structural transition or in the strained state above the 
transition. 

Figure [2] shows images from the strained region B of 
Fig. [I] below (130 K, left top panel) and above (160 K, 
left bottom) the transition. This region is completely de- 
twinned by strain. The red square shows the small clean 
area of the sample, where the color of the image was an- 



alyzed numerically. The right panels in Fig. [2] show red- 
green-blue (RGB) histograms of that area. The images 
were taken every 5 K from 80 to 260 K. Fig. [3] shows the 
difference between the intensities of the blue and green 
channels, indicating that the structural distortion does 
not vanish abruptly at the transition but remains notable 
up to 200 K. 

For reference we show the results of the equivalent 
analysis in the parent compound Cal22. Here the tran- 
sition is strongly first order, and the data show no tail 
above the transition. 

X-ray diffraction 

X-ray analysis was done in both the unstrained (area 
A) and strained (area B) parts of the same crystal, as 
shown in Fig. [I] In the unstrained region, the tetragonal 
(220) Bragg peak splits below Tto into four peaks, each 
representing a distinct orthorhombic domain. The four 
orthorhombic reflections merge back into a single tetrag- 
onal Bragg peak on warming above Tto (see Fig. 4). 

The integrated intensities of the reflections at 6 K allow 
us to determine the relative population of orthorhombic 
domains. In the unstrained area, the population of the 
four domains ranges 19 to 31 % of the total peak inten- 
sity, characteristic of a near random distrobution. In the 
strained region, the domain whose ao-axis lies along the 
direction of the strain accounts for nearly 90 % of the 
probed sample volume. 

These effects are schematically described in Fig. [5] The 
separation between the Bragg peaks resulting from the 
orthorhombic 01 and 02 domains is fixed because the 
relative angle between their twinning planes is fixed. The 
same is true for the peaks from the 03 and 04 domains. 
However there exists no such rule for the separation be- 
tween the 02 and 04 peaks because the angle between 
their twinning planes is determined by their relative do- 
main populations. In a sample with perfectly equal do- 
main populations the four Bragg peaks would produce 
a square with each reflection having equal intensity. As 
the relative population of one domain orientation grows, 
the angle between the twinning planes of the 02 and 04 
peaks becomes smaller and consequently the separation 
of the their Bragg peaks diminishes. This behavior is 
readily seen in the X-ray data in Fig. [4j The unstrained 
region of the crystal manifests relatively similar popula- 
tions of each domain orientation and produces a pattern 
not quite square but slightly trapezoidal below the tran- 
sition temperature. By contrast in the strained region 
of the crystal, where the dominant domain orientation 
represents nearly 90% of the sample volume, the 02 and 
04 peaks are no longer distinguishable as two separate 
reflections. 

The temperature evolution of the orthorhombic dis- 
tortion, e = | a °~^°) ? can be clearly seen as an increased 
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FIG. 6. Orthorhombic distortion, e = \ a ° , , vs. temper- 
ature in strained and unstrained parts of the sample. The 
strain notably increases the orthorhombic distortion below 
the transition, and induces a "tail" of orthorhombic distor- 
tion above the sharp drop at 135 K. 



splitting distance between the orthorhombic reflections 
at 6 K as compared to the splitting at 132 K (see Fig. [4]), 
and is shown in Fig. [6j Application of strain notably in- 
creases e below the transition, and most importantly, a 
"tail" of the orthorhombic distortion can be tracked to at 
least 150 K, well above 135 K where the order parameter, 
e, shows a sharp drop. 



FIG. 7. (Color Online) Temperature-dependent normalized 
resistivity, p a (T)/p a (300 K) of the BaFe2As2 sample in the 
free-standing, p t (black curve), and strain-detwinned, p a (red 
curve) regions of the same sample used for X-ray measure- 
ments in Figs [4] [6] and [9] The third (blue) curve shows pi , 
calculated as pi = 2 * p t — p a • The anisotropy can be seen for 
all temperatures below the transition, and a slight anisotropy 
can be found above the transition. Inset: Progression of the 
effect of increasing strain on the resistivity (p a in the de- 
twinned state). The black curve represents a free standing 
crystal. Tensile stress incrementally increases until reaching 
approximately 20 MPa for strain 5, see text for details. Strain 
2 is sufficient to detwin the sample, revealing a sharp drop in 
resistivity at the transition. On further strain increase the 
jump rounds and its onset shifts up in temperature. 



DISCUSSION 



Resistivity 

Figure [7] shows the normalized temperature-dependent 
resistivity in the twinned, p t: and strain-detwinned, p a , 
states of the same sample as measured by x-ray diffrac- 
tion (Fig [4j [6| . The third curve was calculated assum- 
ing that p t represents an equal mixture of p a and p^, 
Pt =2p t -p a . 

After sufficient stress was applied to detwin the crys- 
tals (the sample in Fig. 7 was nearly completely de- 
twinned after the second strain increment, Strain 2, as 
determined by polarized optical imaging), we performed 
a careful study of the effect of additional stress on the 
resistivity anisotropy. The stress, whose magnitude at 
the highest level is estimated to be in the 20 MPa 
range, increases the onset temperature of the resistiv- 
ity anisotropy. However the most dramatic effects of the 
resistivity change are clearly around the point where the 
strain is sufficient to detwin the crystal. 



Strain-induced Anisotropy 

Figure [8] shows a direct comparison of the temperature- 
dependent degree of the orthorhombic distortion, e(T), 
and of the resistivity anisotropy, pb/pa, m the same sam- 
ple of BaFe2As2 under identical strain conditions. The 
two quantities reveal a clear correlation. Both show a 
rapid rise below approximately 135 K with decreasing 
temperature. In addition both e and pb/pa show a clear 
"tail" above 135 K, in agreement with the color analy- 
sis discussed above. In strain-free samples of Bal22, the 
tetragonal-to-orthorhombic structural transition at Tto 
is of the second order and precedes a strongly first-order 
magnetic transition at T/v < Tto [§]• It is then natu- 
ral to assign the rapid increase of the anisotropy below 
135 K to a magnetic transition, while the "tail" above 
135 K correlates with the orthorhombic distortion. The 
exact meaning of the structural transition in the presence 
of the strain field becomes unclear, as the order param- 
eter varies smoothly with temperature. Therefore from 
this direct comparison we conclude that externally ap- 
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FIG. 8. (Color Online) Comparison of the temperature- 
dependent resistivity anisotropy, y-, and the orthorhombic 
distortion, e = m the temperature range close to 

Tto • Both quantities show a pronounced "tail" above a sharp 
drop in the order parameter at 135 K, revealing that the 
anisotropy is directly related to strain. 



FIG. 9. (Color Online) Comparison of the temperature- 
dependent orthorhombic distortions, e = , in strain 
detwinned areas of SrFe2As2, Ref. [23J and BaFe2As2- The 
data are presented vs. normalized temperature T/Tn- A 
pronounced "tail" above Tn in BaFe2As2 is caused by an 
anomalously strong susceptibility of the lattice to strain. 



plied strain is the cause of the structural and transport 
anisotropy above Tto- 



Comparison of the effect of strain on first and 
second order transition: BaFe 2 As 2 vs SrFe 2 As 2 

In Fig. [9] we compare the temperature dependent or- 
thorhombic distortions, e, for two strained A122 com- 
pounds each with a very different character of the tran- 
sition: strongly first-order in Sri 22 and second order in 
Bal22. The data are plotted on a normalized tempera- 
ture scale, T/Tjsf. As is clear from the Figure |9| the "tail" 
of the anisotropy above the transition is virtually absent 
in Srl22, whereas it is quite noticeable in Bal22. In the 
next section we apply Ginzburg-Landau - type theory 
to model the effect of strain on the resistivity anisotropy 
considering first and second order transitions in the strain 
field. 



Phenomenological model of the effect of the 
uniaxial strain 



Regardless of which electronic degree of freedom cp is 
responsible for the electronic anisotropy, it should be 
proportional to the orthorhombic distortion, since both 
break the tetragonal symmetry of the lattice close to T s . 
For example, cp can be associated with magnetic fluctua- 
tions [32 . By symmetry, and e are bilinearly coupled in 
the free energy expansion, ie. they give rise to the term 
4>e. Since the external strain a also couples bilinearly to 
the orthorhombic distortion e = (ao — ho) /{do + &o)> it 



has an effect on <j) similar to that of a magnetic field h 
on Ising ferromagnets. 

In order to compare the effect of a finite h on the 
second-order and the first-order structural phase tran- 
sitions, we consider the phenomenological free energy: 
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h cp 



(i) 



with temperature parameter r ex T — T s °, where T s ° is the 
mean- field structural transition temperature. Here u and 
w are phenomenological parameters of Ginzburg-Landau 
theory describing the phase transition. To ensure the sta- 
bility of the free-energy expansion, w has to be positive. 
If u > as well, we have a second-order transition. If 
u < 0, we have a first-order transition. In this case, the 
ratio —u/w determines how strong the first-order tran- 
sition is, i.e. what is the magnitude of the jump of the 
order parameter. For u > 0, we have a second-order 
phase transition at r = for h = 0. The effect of a small 
but finite h is to extend the region of finite cp asymptot- 
ically to r — >• 00, giving rise to a "tail" in the plot of cp 



as function of temperature (see Fig. 10). Formally, there 
is no strict T s , although experimentally there will be a 
temperature above which the distortion anisotropy is too 
small to be detected. Notice that, at T s °, the value of (p 
scales with the applied field according to (p ~ h 1 / 6 , where 
S = 3 is the mean- field critical exponent. 

Let us consider u < 0, which gives rise to a first-order 
phase transition. As usual for first-order phase transi- 
tions, there is a coexistence region where the states with 
(p = and (p ^ are both local minima of the free 
energy. If we consider an adiabatic change of temper- 
ature, such that the system always chooses the global 
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FIG. 10. (Color Online) Evolution of the anisotropy parame- 
ter if oc pb — pa vs. temperature parameter r, r oc T — T° . Left 
column of panels is for h = 0, right column is for /i = 0.0 lit;. 
The bottom pair of panels shows a second-order transition for 
u—w—1^ the other pairs of the panels show first-order transi- 
tions for u=-w, u=-0.5w and u=-0.1w (top to bottom). The 
dashed lines show the size of the jump A(p in the absence of 
an external field. 



second order transition originates from the fact that the 
orthorhombic transition is actually not strictly defined 
in the strained (and thus orthorhombically distorted) 
tetragonal phase under uniaxial stress. On the other 
hand it suggests that the susceptibility to stress is no- 
tably enhanced in case of a weak second order transition 
character. 

CONCLUSIONS 

Systematic characterization of the effect of perma- 
nently applied stress on the properties of BaFe2 As2 using 
x-ray, polarized optics and electrical resistivity measure- 
ments suggest that the applied stress is actually the cause 
of the resistivity anisotropy in the nominally tetragonal 
phase. Thus the resistivity anisotropy "tail" above the 
temperature of the structural transition is solely due to 
the effect of the uniaxial strain applied to detwin the sam- 
ples. The difference between AFe2As2 compounds with 
various alkali earth metals is determined by the char- 
acter and the strength (order parameter jump) of the 
structural transition. These conclusions are supported 
by a phenomenological model of the effect of the uniaxial 
strain on the structural transition, similar to the effect of 
a magnetic field on Ising ferromagnets. It is interesting to 
notice that this strong susceptibility to the effect of strain 
is also found in the orbital splitting in ARPES measure- 
ments [HI [223 , and in the temperature dependence of the 
shear modulus in ultrasonic experiments [33] . 
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jump Aip will depend on the value of the ratio h/\u\. 
In Figure [ToJ we plot the temperature evolution of cp for 
different values of \u\ /w — {1, 0.5, 0.1}, keeping h/w = 
0.01 constant. The dashed line shows the magnitude of 
the jump A(f for h = 0. Notice that when the first- 
order transition is stronger, the jump is barely affected 
by the finite field. In particular, above the temperature 
where the jump takes place, the order parameter is never 
zero but is always very small, giving rise to a rather small 
"tail" . On the other hand, when the first-order transition 
is weaker, the same field can completely smear out the 
jump. This gives rise to a noticeable and continuous 
"tail" , and therefore to a second-order transition [34] . 

This analysis suggests that the anisotropy above the 
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